Chiral transition in a magnetic field and at finite baryon density 
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We consider the quark-meson model with two quark flavors in a constant external magnetic field B 
at finite temperature T and finite baryon chemical potential hb- We calculate the full renormalized 
efTective potential to one-loop order in perturbation theory. We study the system in the large- 
Nc limit, where we treat the bosonic modes at tree level. It is shown that the system exhibits 
dynamical chiral symmetry breaking, i. e. that an arbitrarily weak magnetic field breaks chiral 
symmetry dynamically, in agreement with earlier calculations using the NJL model. We study the 
influence on the phase transition of the fermionic vacuum fluctuations. For strong magnetic fields, 
\qB\ ~ 5m^ and in the chiral limit, the transition is first order in the entire )j,b —T plane if vacuum 
fluctuations are not included and second order if they are included. At the physical point, the 
transition is a crossover for /ib = with and without vacuum fluctuations. 

PACS numbers: 11.15Bt, 04.25.Nx, ll.lOWx, 12.38Mh 



I. INTRODUCTION 

A tremendous amount of work has been done in re- 
cent years to map out the QCD phase diagram as a func- 
tion of temperature T and baryon chemical potential fi b ■ 
Two aspects have received particular attention: the posi- 
tion of the critical endpoint where the curve of first-order 
chiral transitions terminates in a second-order transi- 
tion [l[ , and the various color-superconducting phases at 
large baryon chemical potential and low temperature [2] . 

With Nf quark flavors, the global symmetry of QCD 
is SU{Nf)v X SU{Nf)A in the chiral limit and SU{Nf)v 
if the quark masses are equal. In the chiral limit with 
Nf = 2, we use the isomorphism between the groups 
SU{2) X SU{2) and 0(4) and use the 0(4) linear sigma 
model (LSM) as a low-energy effective theory for QCD. 
Coupling the linear sigma model to quarks we obtain the 
quark-meson (QM) model. Since we have quark degrees 
of freedom, we can couple the model to a baryon chem- 
ical potential fiB and study finite-density effects. The 
QM model has been used to study various aspects of 
the chiral transition at /is = j^-|7[ and (8l-[l0|. 

Schwinger-Dyson equations were used in [12l |. One can 
couple the Polyakov loop to the quark sector by intro- 
ducing a constant temporal gauge field background fisl] , 
in order to incorporate certain aspects of the deconfine- 
ment transition. The thermodynamics of the Polyakov- 
extended quark-meson model (PQM) was studied in 
Refs. [H-ia, and the PNJL model in Refs. [H,!!!-!!] 

Using the QM model, it is a common approximation 
to omit the quantum and thermal fiuctuations of the 



* Electronic address: 
t Electronic address: 



andersen@tf . phys . no | 



rashid. khan@ntnu.nol 



bosonic degrees of freedom, i.e. treating them at tree 
level. This is the large-TVc limit. On the other hand, 
one keeps the thermal fluctuations of the quarks. Since 
symmetry breaking takes place in the mesonic sector and 
at tree level, the vacuum fiuctuations of the quarks are 
sometimes omitted; however one must be careful when 
making such additional approximations. It was shown in 
Ref. 7] that the order of the phase transition depends on 
whether or not one includes the fermionic vacuum con- 
tributions in effective potential. This issue was studied 
in detail in [10] using optimized perturbation theory (ll| . 
If all bosonic and fermionic vacuum and thermal contri- 
butions are included, resummed one-loop results predict 
a first-order transition in the entire /i-T plane in the 
chiral limit. 

For Nf — 2 and in the chiral limit, universality and 
renormalization arguments firmly establish the second- 
order nature of the transition at zero baryon chemical 
potential if the axial U{1) symmetry is explicitly bro- 
ken. Otherwise it is driven first order. For Nf = 3, 
the transition is always first order in the limit of zero 
current quark masses. [H, [s^l- At the physical point, i. 
e. for current quark masses that yield — 140 MeV, 
it is a crossover transition. At nonzero baryon chemi- 
cal potential and in particular at T = 0, the order 
of the chiral transition is not obvious from universality 
arguments [3l[. Most model calculations predict a first- 
order transition at T = ^, ^] , although the results 
depend on the input parameters and are sensitive to the 
ultraviolet cutoff. If there is a first-order transition at 
r = 0, there is a line of first order transitions in the 
/is-T plane that ends at critical point since we known 
that the transition is in a second order for ^b = 0. 

Another important question is the behavior of QCD 
in external magnetic fields. This is not a purely aca- 
demic question since the properties of QCD in strong 



magnetic fields B is relevant in several situations. For 
example, large magnetic fields exist inside ordinary neu- 
tron stars as well as magnetars [s^. In the latter case, 
the cores may be color superconducting and so it is im- 
portant to study the effects of external magnetic fields in 
this phase. Similarly, it has been suggested that strong 
magnetic fields are created in heavy-ion collisions at the 
Relativistic Heavy-Ion Collider (RHIC) and the Large 
Hadron Collider (LHC) and that these play an impor- 
tant role [S^l- In this case, the magnetic field strength 
has been estimated to be up to i? ~ 10^^ Gauss, which 
corresponds to \qB\ ~ 6m^, where \q\\s the charge of the 
pion. One of the interesting issues that has been raised 
is whether there is a splitting of the deconfinement and 
chiral transitions at finite B. This has spurred the inter- 
est in studying QCD in external fields. At zero baryon 
chemical potential this can be done from first principles 
using lattice simulations [ssl - ISTj . At finite this is very 
difficult due to the infamous sign problem. Therefore 
one often resorts to effective theories that share some of 
the features of QCD, such as chiral symmetry breaking. 



A low-energy effective theory that provides a sys- 
tematic framework for systematic calculations is chi- 
ral perturbation theory. Chiral perturbation theory 
has been used to study the quark condensate in strong 
magnetic fields at zero temperature [H, H^, while 
the quark-hadron phase transition was investigated in 
Ref. 40] . The effects of external magnetic fields on the 
chiral transition have been studied in detail using the 
NJL model klMSOl . the Polyakov-loo p e xtended NJL 



model 
model 



53, the QM model [49|, |50|, |5J] , the (P)QM 
and the instanton-liquid model [s^. Simi- 
larly, the effects of strong magnetic field on the various 
color superconducting phases have been studied using 
the NJL model ^s-M- 



In this paper, we study the effects of magnetic fields 
on the chiral phase transition at finite temperature and 
finite baryon density. In order to do so, we employ the 
quark-meson model with two fiavors coupled to a baryon 
chemical potential /is. We consider the large iVc-limit 
which amounts to treat the sigma and pious at tree level. 
In this limit, we investigate the role of the fermionic 
vacuum fluctuations. 



The paper is organized as follows. In Sec. II, we briefly 
discuss the QM model in an external magnetic field and 
at finite ^xb- In Sec. Ill, we calculate the standard one- 
loop effective potential. In Sec. IV we discuss the ther- 
modynamics and the numerical results in the large- iVc 
limit. In Sec. V, we summarize and conclude. 



II. QUARK-MESON MODEL 

The Euclidean Lagrangian for the quark-meson model 
with Nf — 2 fiavors is 

^ — -^mcson -^quark -^Yukawa ~^ ^dct ; (1) 

where the various terms are 

^meson = Tr [9^$^^^$] + m^Tr 

-f-Tr[$t$]^_i/iTr[$ + $t] , (2) 
3 2 



>C Yukawa = 94' W - «75T • Tt] Ip 

£dot = cdet[$-l-$t] ^ 

where 

1 



{(7 + T ■ Tv) 



(3) 
(4) 
(5) 

(6) 



Here, a is the sigma field, tt denotes the neutral and 
charged pions. Moreover t are the Pauli matrices, fi = 
^(Mji + Mrf) is the quark chemical potential, where /i„ 
and fid are the chemical potential for the u and d quarks, 
respectively. The baryon chemical potential is given by 
fiB = 3/i. We set fiu = fJ-d so that we are working at 
zero isospin chemical potential, /i/ — ^(/iu — /id) = 0. 



The Euclidean 7 matrices are given by 



«7m 



and 



74 = 7^^^, where the index M denotes Minkowski space. 
The fermion field is an isospin doublet 



4 



(7) 



li h = 0, the first three terms in Eq. ([T]) are invariant un- 
der U{2)l x 1/(2)^ ~ SU{2)l x SUi2)R x ;7(l)s x 
If /i ^ 0, chiral symmetry is explicitly broken, otherwise 
it is spontaneously broken in the vacuum. Either way, 
the symmetry is reduced to SU{2)v x U{1)b x U{1)a- 
Note that this requires < which is assumed in the 
remainder of the paper. The U{1)a symmetry is also 
broken in the vacuum by instantons |62| . and their ef- 
fects are mimicked by the the determinant term Eq. ([5]) . 
Since the U{1)b symmetry is always respected, symme- 
try of the QCD vacuum is SU{2)v x U{1)b- In the 
following we set c = for simplicity. 

Chiral symmetry is broken in the vacuum by a nonzero 
expectation value v for the sigma field. We therefore 
make the replacement 

a V + a , (8) 

where a is a quantum fiuctuating field with vanishing 
expecation value. After the shift the tree- level po- 
tential is given by 



^v^-hv + ^B'+£o, (9) 
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where the last term is the vacuum energy density . The 
tree-level masses for the sigma and neutral pion are 



2 I 2 

m + —V 

2,-^2 

m H V 

6 



(10) 
(11) 



Note that the pion mass vanishes at the minimum of 
the tree-level potential in accordance with Goldstone's 
theorem. The charged pions tt^ and the quarks couple 
to the external magnetic field. This coupling is imple- 
mented by the substitution ^ + iq^fi, where 
is the four- vector potential and q is the electric charge 
of the particle. For a constant magnetic field B in the 
z-direction, one can conveniently choose the four-vector 
potential as {Aq,A) — {0,0, Bx,0). The classical solu- 
tions to the Klein-Gordon equation in a constant mag- 
netic field are well known and the dispersion relation is 
given by 

iKpf = pI+ + \^^^ + (2" + 1)19^1 - (12) 

where n is an nonnegative integer, q is the electric charge 
of the pion, and Pz is the spatial momentum in the 
z-direction. The subscript ± denotes tt^ and we note 
that the dispersion relations are identical. Similarly, the 
Dirac equation in a constant magnetic field B can be 
solved straightforwardly and the dispersion relation for 
the quarks is given by 



El 



Pz 



+ {2n+l^ s)\qfB\ 



(13) 



where niq — gv is the quark mass after symmetry break- 
ing, g/ is the electric charge of the quark, and s = ±1 
denote spin up/down, respectively. 



III. ONE-LOOP EFFECTIVE POTENTIAL 

In this section, we calculate the one-loop effective po- 
tential of the quark-meson model. This is done by taking 
into account the Gaussian fluctuations around the mean 
field V. The one- loop contribution Vi to the effective 
potential can be written as the sum of the contributions 
from the cr, pions, and quarks. This yields 



Vl = V, + V^O + + V^- + Vq , 



where 



V. = ^^log[P'+m^] 
= log -f m^] 



(14) 

(15) 
(16) 



= / log[Po^+P^ + AfI] , (17) 



f 



rrio 



(18) 



where = m'^ + {2n + l)\qB\ and / denotes the flavor. 
The symbol is short-hand notation for 



in 



T 



Pa=2-nnT 



(19) 



where d — 3 — 2e and A is the renormalization scale 
associated with dimensional regularization in the MS 
scheme. The integral over spatial momenta will be cal- 
culated in dimensional regularization. In the case of the 
charged pions, the sum-integral is replaced by a sum over 
Matsubara frequencies Pq = ^irnT, a sum over Landau 
levels n, and an integral over momenta in d — 2 = 1 — 2e 
dimensions: 



\qB\T (e^^K 



2-K 



An 



E 



(20) 



where the prefactor takes into account the degen- 
eracy of the Landau levels and 



{2nY- 



(21) 



Similarly, the contribution from the quarks can be writ- 
ten as 



E 

{Po},n,s,f 



\qfB\ 

2n 



log [PS 



(22) 



where = 



r2 _ _2 , (2n + 1 - s)\qfB\, Nc is the number 
of colors, and the Matsubara frequencies are Pq = {2n + 
l)nT. 

The sum-integrals involving the contribution from the 
sigma and the neutral pion are essentially the same so 
we consider Eq. (|15p with a general mass M. Summing 
over Matsubara frequencies we can write 



^log [P' 



M 



/r21 



+2Tlog 



■ (23) 



The first term in Eq. is ultraviolet divergent. Calcu- 
lating it with dimensional regularization and expanding 



3 



in powers of e through order e*^, we obtain 



i 



log 



T 





f-V 


'1 


3" 


327r2 


{m^ J 


e 


^2 



dpp^ log 



1 



.(24) 



We next consider the contribution to the effective po- 
tential from the charged pions. Summing over the Mat- 
subara frequencies in Eq. (jl8p . we obtain 



47r 



E 



Pi 



-2rioK 1-e 



(25) 



The first integral is ultraviolet divergent and we compute 
in dimensional regularization with d = 1— 2e. This yields 



Ml 



Ml 
47r 



r(-l + e) (26) 



Eq. ([261) shows that the sum over Landau levels n in- 
volves the term M^~^'. This sum is divergent for e = 
and we regulate it using zeta-function regularization. 
After scaling out a factor of {2\qB\Y^'^, this sum can 
be written as 



^M^^ = {2\qB\f-^Y: 



2|gB|_ 

{2\qB\f-%{-l + e,\ + x) , (27) 



where x = ^™^| and ('(q, s) is the Hurwitz zeta function. 
The vacuum contribution then reduces to 



Vvac 



87r2 



im) r(-i + ^K(-i + e,H-)- 

(28) 



Expanding Eq. (|28|) . we obtain 



^ 1 f A2 yifjqBy 

647r2 \2\qB\J [\ 3 



+ 1 + 8(gB)2c(i'")(-l, i + a;) + 0(e) 



(29) 



where C^^'°H-1: ^ 

-|- x) is the derivative of the Hurwitz zeta function with respect to the first argument and where 

we have used that C(— 1: | + ^) — ^ 

The vacuum contributions from the quarks can be calculated in the same manner, and one finds [46l [63l [63 | 



-\;vac _ c 
« " 2^2 



2\qiB\ 



l&n^^\2\qfB\) 



C(-l + e,a;/)- ^4 
f 



- + 1 ) - 8{qfByC''''^i-l,Xf) - 2\qfB\m'^\ogXf + 0{e) 



(30) 



where Xf = m'^/2\qfB\. The divergences of the effective 
potential are given by Eqs. and The 

divergences that depend on the magnetic field are given 
by 



Eq. ® [H: 



B' 



B' 



1 - 



487r2e 



f27r2e 



(32) 



[qBf 
967r2e 



247r2e 



(31) 



These divergences are removed by wavefunction renor- 
malization of the (external) gauge field A^. This is done 
by making the replacement in the tree-level Lagrangian 



The remaining divergences in Vi are given by 
1 , 



idiv 



647r2e 



These are the same divergences as one encounters in van- 
ishing magnetic field and so the usual renormalization 
procedure can be used to eliminate them. This is done 



4 



by adding a vacuum energy counterterm A£q and mak- _ _ •iNcNfg'^ . . 

ing the replacements m? + 5m? and A ^ A + 5 A Syr^e 27r2e 

in the tree-level effective potential where 

= ' (33) 

The renormalized one-loop contribution to the effective 
5m?' = — , (34) potential at T = then becomes 



log 



A2 



647r2 

_M)!c(^^°)(-i - 

(gi?)2 A2 
log ■ 



647r2 



log 



A2 



327r2 



log 



A2 
2|gi?| 



167r2 



E 



log 



A2 
2|9/S| 



967r2 2|gB| 



2\qfB\ 



log 



A2 



247r2 2\qfB\ ' 



(36) 



The finite-temperature contribution v'[''^ to the effective potential is given by 



T 



dpp log 



1-e 



T 



dpp log 



1 - e 



|gS|T. 



E / '^p 



1 - e 



\lfB\T 
27r2 



dp log 



(37) 



The full one-loop effective potential is then given by the 
sum of Eqs. ([9]), (|36|) . and (l37l) . The renormalized vac- 
uum energy background field B, mass m, and quartic 
coupling A satisfy the renormalization group equations 



A 



A 



'dA 
dB^ 
'dA 

dm^ 



m X 
8^ ' 

247r2 



67r2 



A 



dA 
dA 



8^2 ' 

A2 mcNfg^ 



dA 



167r2 



(38) 
(39) 

(40) 
(41) 



Using the renormalization-group equations, we conclude 
that the one-loop effective potential is renormalization 
group invariant. 

We close this section by taking the limit B in 
Eqs. (|36)) and (|37)) . We then need the large- a; behavior 
of ((1.0) (-1,1 +2;) and C^i'")(-l,a;/). These are given 



by 



((1.0)(_i 



log a: 



C(1'0)(-1,.t;) = -x} [\ogxj-- 




(42) 



-xf \ogXf + 



(43) 



Inserting the expansions (|42|) and (143)) into ([361) . we ob- 
tain the standard one-loop vacuum term 



647r2 



r A2 


3' 


3mi 




A2 


3" 


log — - 






log 








^2 


647r2 




^2 



167r2 



, A2 3' 
log + o 



(44) 



In Eq. (P7)) . we change variable = 2|gi?|n, which 
yields p±dp± = \qB\dn. Replacing the sum by integrals, 
we obtain 



5 



n^'" = / dp log 



1 - e 



ST r ^ 2, 

2^ X "^"^ 



1 _ e-^\/p^ 



dp p log 



1 j_ g-;3(^p2+,„2±^,) 



(45) 



The full one-loop effective potential at i? = is the sum of Eqs. Q, and (|^. 



IV. THERMODYNAMICS AND NUMERICAL 
RESULTS 

In the previous section, we calculated the one-loop 
effective potential at finite temperature T and finite 
baryon chemical potential /is. A common approxima- 
tion in the QM model is to neglect the quantum and 



thermal fluctuations of the mesons, which is equivalent 
to the large- iVc limit. One hopes that the contribu- 
tions from the quarks include the most important ef- 
fects [3, '53]. In this paper, we will apply this approxi- 
mation and defer the use of more sophisticated methods 
to a subsequent paper [g^. The one-loop effective po- 
tential then reduces to 



Vo-i-i = -m^v^ + —v^ ~ hv 



Ncml 
87r2 



E log 



NcTTl 

^ 167r2 
2\qfB\ 



/ 



log 



A2 



2\qfB\ 

kfBjT 
27r2 



N, 



dp log 



(46) 



where we for simplicity have omitted all ^-independent 
terms. Notice that the renormalized effective poten- 
tial Eq. for large values of the order param- 
eter V is unbounded from below due to the domi- 
nant term —Ncm^\og{m^/2qfB)/16TT^ that arises from 
^(i,o)(_1^2;^) at large Xf. The unboundedness of the 
fermionic functional determinant is well known for B = 
and arises here as well. It is due to the minus sign 
in Eq. ([T8|) . However, it is stabilized by the bosonic 
contributions in Eq. (|44l) for realistic values of the cou- 
pling A, (cf. Eqs. (|49l) - ([52l) and the numerical values 
given below). Generally, the reliability of perturbative 
(vacuum) calculations and the stability of the effective 
potential are nontrivial issues, see e. g. for a thorough 
discussion I67H69I. 



field B corresponds to choosing > ^. The mini- 
mum of the effective potential is found by minimizing 
the effective potential with fj,B = T = h = 0. In 
the limit S — > 0, we know that the expectation value v 
is exponentially small [6^ and so we can expand Vo+i 
around v = 0. This yields ^ 



dV, 



0+1 



dv 



E \if^\ 



Either u = or we obtain v'^ 



2i\qB\ 



exp 



Nc9^\qB\ 



(47) 



(48) 



The nontrivial solution Eq. (|48)) corresponds to the local 
minimum of Vq+i- Eq. (1481) has the same functional form 



Let us first briefly discuss magnetic catalysis of dy- 
namical symmetry breaking, namely the effect that chi- 
ral symmetry is broken dynamically for any nonzero 
magnetic field when it is intact for B = 0. This ef- 
fect is now well established, see e. g. [gl, [70l - [75| . In 
the present case, chiral symmetry at vanishing magnetic 



^ > corresponds to a symmetric vacuum state at tree level. 
In the NJL model, symmetry breaking is always a loop effect 
which takes place for G smaller than a critical coupling Gc- 

^ We have used (^^^''^^{—l,x-f) = const — x f log\/2TTX f + ... for 
small values oi Xf. 



6 



with the replacement l/G — > rv? j as in the NJL-model 
calculations of e. g. Ref. [63]. The reason is that we 
are basically evaluating the same fermionic functional 
determinant although we are using a different ultraviolet 
regulator A. 

We next discuss the determination of the parameters 
of the Lagrangian ([1]). In the vacuum, we have fiB = 
and /jr = where is the pion decay constant. At tree 
level, the mass parameter m^, the couplings A and g, and 
the symmetry breaking parameter h can be expressed in 
terms of observable masses of the sigma and pions, and 
the pion decay constant fj^: 



m? — 


1/ 2 
{mt 


A = 








ruq 


.9 = 


U ' 


h = 





^2 \ 



(49) 

(50) 

(51) 
(52) 



We use a sigma mass of = 800 MeV, a constituent 
quark mass of mq = 300 MeV, and a pion decay constant 
of /tt = 93 MeV. In the chiral limit, h = 0, this yields 
to2 = -320000 MeV^, A = 222, and g = 3.2258. At 
the physical point, we obtain to^ — —291018 MeV^, A = 
215.29, and g = 3.2258. In the remainder of the paper, 
we set A^c = 3 and c = 0. 

We first consider the effective potential at T = /j, = 
0. As noted above, we have determined the parameters 
of the theory at tree level. Quantum fluctuations, i. 
e. loop corrections, will modify the classical potential 
and it depends on the renormalization scale A. We can 
choose A such that the minimum of the one- loop effective 
potential in the vacuum (for B = 0) still is at v = Jt^. 
This is done by requiring 



0+1 



= . 



(53) 



This is equivalent to requiring that the one-loop self- 
energy of the pion at zero external momentum vanish, 
ni,^(0) = 0. From Eq. dM]), this yields 



log 



A^ 



1 



= 0, 



(54) 



whose solution is A = 181.96 MeV. We will use this value 
in the remainder of the paper. 

In Fig. [U we plot the tree-level potential (solid curve) 
as well as the one-loop effective potential in the chiral 
limit for \qB\ = 0, \qB\ = 5mJ, and \qB\ = lOm^ with 
T — fjL = {) The dashed curve is for \qB\ = 0, the dot- 
ted curve is \qB\ = 5toJ and the dash-dotted curve is 
for \qB\ = lOm^. We notice that the potential becomes 



deeper with increasing magnetic field and that the min- 
imum moves to larger values oi v ^ . Thus the magnetic 
field enhances symmetry breakin g. T his is in agreement 
with earlier findings [Hj, [H, l63l, l70l - l75j . We also note 




V [MeV] 

FIG. 1: Normalized effective potential Vo+i/fn in the chrial 
limit for T — hb ~ 0. Tree level (solid curve), one-loop with 
\qB\ = 0(dashed curve), one loop with \qB\ — 5m^ (dotted 
curve), and one loop with |g_B| = lOm^ (dash-dotted curve). 

that the local maximum of the effective potential for 
large values of v becomes lower with increasing mag- 
netic field. Thus for sufhciently strong magnetic fields, 
the local minimum at = 93 MeV ceases to exist and 
the system has no longer a metastable state. This is the 
instability that we discussed above. This in disagree- 
ment with Ref. [55'[ . In their work, the authors subtract 
the fermionic vacuum fluctuations at -B = 0. Thus they 
subtract a w-dependent contribution given by the last 
term in Eq. (|44l) which leads to a one-loop effective po- 
tential that is independent of the renormalization scale 
A and is stable for large values of v. 

It was mentioned in the introduction that the inclu- 
sion of the fermion vacuum fluctuation term can change 
the order of the phase transition in the chiral limit as 
well as strongly influence physical observable. This was 
discussed in detail for = in Ref. Q- It turns out 
that this is the case in the entire T-^b phase diagram 
for B = [l3[. We find that this is also the case in a 
strong magnetic field. In Fig. [21 we show the normalized 
effective potential with no vacuum fluctuations for three 
different values of the baryon chemical potential /is and 



Note, however, that the one-loop effective potential for \qB\ = 
is more shallow than the tree-level potential. 
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T = as a function of v. Clearly, the transition is first 
order if the vacuum fluctuations are not included. 



10 - 



5 - 
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20 40 60 80 100 120 

V [MeV] 

FIG. 2: Normalized effective potential Vo+i/f^ without 
vacuum fluctuations for T = and \qB\ = 5m^ as a function 
of V for three different values of /x: fi — (dashed curve), 
/i = /ic = 301.6 MeV (sohd curve), and = 358 MeV (dotted 
curve) . 

Similarly, in Fig. [31 we show the normalized effective 
potential with vacuum fluctuations for three different 
values of the baryon chemical potential /is and T = 
as a function of v. The transition is second order if the 
vacuum fluctuations are included. 

Vo+i//^ 



15 - 




-iq L 7 r J 

20 40 60 80 100 120 140 

V [MeV] 

FIG. 3: Normalized effective potential Vo+i/f* including 
vacuum fluctuations for T = and \qB\ = 5m^ as a function 
of V for three different values of /i: fi — (dashed curve), 
fj, = fj.^ = 324 MeV (solid curve), and fi = 380 MeV (dotted 
curve) . 

In Figs.HHil we show the phase diagram as a function 
of fiB and T in the chiral limit. In Eq. pS)) . there is 



a sum over Landau levels and for each value of /^b, we 
include sufficiently many terms so that our result for 
Tc is converging. For example, for fis = 0, we must 
typically sum the first ten term in the series, while for 
T = it suffices to include the first few terms in the 
series. 

In Fig|4l we have included the fermionic vacuum fiuc- 
tuations and the dashed curve indicates a second-order 
transition. In the case where we include the quan- 
tum fiuctuations, the critical temperature for /i = is 
Tc — 179 MeV and the critical chemical potential for 
T = is /ic = 324 MeV. In comparison, the critical tem- 
perature at vanishing magnetic field is Tc = 176 MeV, 
the critical chemical potential is /ic = 321 MeV, and 
the phase transition is also second order for all values of 
H [l3]- Calculating the critical temperature fj,B = 
and \qB\ = 10m J yields Tc = 184 MeV and so it seems 
to be increasing very weakly with the strength of the 
magnetic field. 

In Fig [5l we have omitted the vacuum term and the 
solid curve indicates a first-order transition. In the case 
where the quantum fiuctuations are omitted, the critical 
temperature for /i = is Tc = 160 MeV and the critical 
chemical potential for T = is /ic = 301.6 MeV. In com- 
parison, the critical temperature at vanishing magnetic 
field is Tc = 179 MeV, the critical chemical potential 
is /ic — 347 MeV, and the phase transition is also first 
order for all values of /i [l(l|. In contrast to the zevo-B 
case, the exclusion of the vacuum fiuctuation decreases 
the critical temperature significantly. 

In Ref. |10|[. it was shown using optimized perturba- 
tion theory |ll| that the inclusion of the vacuum and 
thermal contributions from the bosons significantly low- 
ers the temperature of the transition. Whether this im- 
plies that including the effects of the bosons will lower 
the critical temperature for nonzero B as well is not 
known. In particular, the i?-dependence of the critical 
temperature is an open question, see also the discussion 
below. 

There have been a number of model calculations us- 
ing chiral perturbation theory, the NJL model, the NJL 
model coupled to the Polyakov loop, the QM model, 
and the Polyakov-loop extended model. Except for chi- 
ral perturbation theory [i^l , they all predict an increase 
of the critical temperature with increasing magnetic 
field. This is in disagreement with the lattice simula- 
tions of Ref. [13] which show a significant decrease of 
the crossover temperature. 

In Ref. (ioj . the authors are using chiral perturbation 
theory to investigate the quark-hadron phase transition 
as a function of the magnetic field at the physical point. 
They compare the pressure in the hadronic phase with 
that of the quark-gluon plasma phase for an ideal gas 
of quarks and gluons, and subtracting the vacuum en- 
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T [MeV] 



HB [MeV] 

FIG. 4: Phase diagram in the fiB-T plane for — 5m^ in 
the chiral limit. Fermionic vacuum fluctuations are included. 

r [MeV] 



fiB [MeV] 

FIG. 5: Phase diagram in the fiB-T plane for \qB\ = 5m^ in 
the chiral limit. Fermionic vacuum fluctuations are excluded. 



ergy due to a nonzero gluon condensate {g'^G^^G^") . 
For weak magnetic fields, the transition is first order. 
Tlie line of first-order transitions ends at critical point 
iV\qB\,T) = (600, 104) MeV. For larger values ofja-Bj, 
the transition is a crossover. The authors of Ref. '55| re- 
port that they find a first-order transition at the physical 
point if they keep the remaining i?-dependent vacuum 
fluctuations and a crossover if they are ignored. Our 
sample calculations at the physical point for \qB\ = brn^ 
(not shown) suggest that the transition is a crossover (for 
Mb = 0). This is in agreement with the lattice simula- 
tions in Ref. [s^j. The disagreement between Ref. [ssj 
and the present results can probably be traced back to 



a different treatment of the vacuum fluctuations, as dis- 
cussed earlier. The disagreement between various ap- 
proaches deserves further investigation. 

V. SUMMARY AND OUTLOOK 

In the present work, we have calculated the one-loop 
effective potential for the quark-meson model at finite 
temperature and baryon density in an external magnetic 
field. We have made the common approximation where 
we ignore all quantum and thermal effects of the bosons 
and hence treat them at tree level. We have seen that 
the critical temperature and the order of the phase tran- 
sition depend on whether one includes the vacuum fluc- 
tuations. Taking a model seriously means including the 
effects of all its degrees of freedom. It therefore ques- 
tionable to throw out some terms unless one can show 
that they are not important. For example, the effective 
potential is receiving contributions from vacuum fluctu- 
ations from all energies scales up to the ultraviolet cutoff 
A of the effective theory. This is the region of validity 
of the low-energy effective theory and it is therefore not 
obvious that they can be omitted. 

Given the conflicting results of various model calcula- 
tions and lattice calculations, clearly more work needs 
to be done. On the analytical side, an important exten- 
sion is to include the vacuum and thermal effects of the 
bosons. The dominant contributions at weak coupling 
from the bosons arise from the daisy or ring diagrams. 
There has been some work on resumming the ring di- 
agrams in the presence of a magnetic field. For exam- 
ple, the authors of Ref. [zlj investigated their role in 
the electroweak phase transition in the standard model. 
Similarly, their effects at weak fields have been studied 
in Ref. [73] in the context of the chiral transition. One 
should apply sophisticated resummation techniques or 
nonperturbative methods such as optimized perturba- 
tion theory ill], the 2PI effective action formalism 781 ] . 
or the functional renormalization group 79]. These 
methods are known to correctly predict a second-order 
phase transition for /is = = [§, [13, HH • 
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